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. The quantum-mechanical collapse (alias fall onto the center of particles attracted by potential 

^s^j ■ ~ r ~ 2 ) is a well-known issue in the quantum theory. It is closely related to the quantum anomaly, i.e., 

breaking of the scaling invariance of the respective Hamiltonian by the quantization. We demonstrate 
, that the mean-field repulsive nonlinearity prevents the collapse and thus puts forward a solution 

• to the quantum-anomaly problem different from that previously developed in the framework of 

the linear quantum-field theory. This solution may be realized in the 3D or 2D gas of dipolar 
bosons attracted by a central charge, and in the 2D gas of magnetic dipoles attracted by a current 
filament. In the 3D setting, the dipole-dipole interactions are also taken into regard, in the mean- 
field approximation, resulting in a redifinition of the scattering length which accounts for the contact 

^ , repulsion between the bosons. In lieu of the collapse, the cubic nonlinearity creates a 3D ground 

state (GS), which does not exist in the respective linear Schrodinger equation. The addition of the 

bJT)' harmonic trap gives rise to a tristability, in the case when the Schrodinger equation still does not 

1 1 1 | lead to the collapse. In the 2D setting, the cubic nonlinearity is not strong enough to prevent the 

collapse; however, the quintic term does it, creating the GS, as well as its counterparts carrying 
the angular momentum (vorticity). Counter-intuitively, such self-trapped 2D modes exist even in 
the case of a weakly repulsive potential r~ 2 . The 2D vortical modes avoid the phase singularity at 
the pivot (r = 0) by having the amplitude diverging at r — > 0, instead of the usual situation with 

+_> ■ the amplitude of the vortical mode vanishing at r — > (the norm of the mode converges despite of 

the singularity of the amplitude at r — > 0). In the presence of the harmonic trap, the 2D quintic 
model with a weakly repulsive central potential r -2 gives rise to three confined modes, the middle 
one being unstable, spontaneously developing into a breather. In both the 3D and 2D cases, the GS 
wave functions are found in a numerical form, and also in the form of an analytical approximation, 

C , which is asymptotically exact in the limit of the large norm. 

, PACS numbers: 03.75.-b; 03.75.Kk; 05.45.Yv; 03.75.Lm 
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> ' I. INTRODUCTION AND THE PHYSICAL SETTINGS 

o : 

It is well known that attractive potential U(r) — — (Uo/2)r~ 2 plays a critical role in quantum mechanics. Indeed, 
(/"■) ■ while the corresponding Hamiltonian, taken in the scaled form, H = — (1/2) V 2 + U (r), obeys the scaling invariance, 

H(ar)=H(r)/a 2 , (1) 

o : 

, the quantization breaks the invariance, which is known as the quantum anomaly , alias "dimensional transmutation" 
A manifestation of the anomaly is that, at Uq > (Uq) 1 ^^ = 1/4, the corresponding 3D (three-dimensional) 
Schrodinger equation, iip t = Hip, does not produce the ground state (GS); instead, it gives rise to the quantum 
collapse, also known as fall onto the center Q. In the 2D space, the critical value of the potential's strength is 

exactly zero, (Uo) cl . = 0, and in the ID setting the same potential gives rise to a still stronger supers election effect, 
effectively splitting the ID space into two disjoint subspaces, i^O ||. 

A solution to the quantum-anomaly problem in the 3D case was proposed outside of the realm of quantum mechanics, 
in terms of a linear quantum-field-theory description, which introduces a renormalization and imposes a GS with an 
arbitrary spatial scale [l|. The purpose of the present work is to demonstrate that a different solution to this 
problem can be obtained within the framework of the mean-field approach, taking into account the effective repulsive 
nonlinearity induced by collisions of particles trapped in the potential under the consideration. In the 3D setting, the 
usual cubic self-repulsive term will be sufficient for this, while in the 2D geometry a stronger quintic term is necessary. 
As a result, we will find the GS with a scale uniquely defined by the physical parameters. It will be demonstrated that 
the same nonlinearity is helpful too below the critical point, where it makes the GS normalizable [at Uq < {Uo) CT ' , 
the linear Schrodinger equation gives rise to GS wave functions with the divergent norm, see below] . 
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The linear Schrodinger equation with the critical potential is 

Wt = -\ (V 2 + U r- 2 - flV) V, (2) 

where the external trapping potential, (J7 2 /2) r 2 , is included too [4]. A physical realization of Eq. in the 3D space 
is provided by molecules with a permanent electric dipole moment, d, interacting with charge Q placed at the origin, 
which creates electric field E = Qr/r 3 . As demonstrated in a recent experiment, the charged particle (ion) immersed 
into an ultracold gas may be kept at the central position by means of the laser-trapping technique @ • Assuming that 
the orientation of the dipole is locked to the local field, i.e., d = sgn(Q)d(r/r), the respective interaction potential is 
U(r) = d ■ E, which corresponds to Eq. ([2]) with Uq — 2\Q\d. This realization is relevant to the capture of electrons 
by dipolar molecules, which was studied in detail experimentally [f| 0], and to Bose-Einstein condensates (BECs) 
formed by dipolar molecules, such as Li-Cs We stress that we consider the gas of permanent dipoles, while, in 
the case of polarizable molecules, with d = eE, where e is the polarizability, the attractive potential is different, 
U = -eQ 2 r~ 4 . 

If the gas of ultracold dipolar molecules is trapped in a pancake-shaped configuration sustained by an appropriate 
external potential [9(, inserting the central electric charge provides for the realization of the 2D version of Eq. ©. 
An additional realization of the 2D setting is offered by a gas of atoms (chromium [l(| [ll|) or molecules (such as 
87 Rb2 [l3j) carrying magnetic moments, the attractive potential being induced by the magnetic field of a transverse 
current filament, with the orientation of the dipoles locked to the local magnetic field. Actually, the orientations of 
the dipoles may form the configurations obeying Eq. @ in the dipolar BEC produced by means of the all-optical 
trapping [11], which does not freeze the dipoles into the confining magnetic field. 

In the mean-field approximation, the contact repulsive interaction in the bosonic gas is represented by the cubic 
term [12]. With the addition of this term, the linear Schrodinger equation ^ is replaced by the Gross-Pitaevskii 
equation (GPE), 

ii/H = - (1/2) (V 2 + U r- 2 - fiV) t/, + i>- (3) 

The relation between the scaled variables and constants, in terms of which Eq. ([3]) is written, and their counterparts 
defined in the usual physical units [12| is 

r = — > 1 = — *ph. V- = 2yi^7roV P h, U = (U Q ) ph , Q = -^Op hj (4) 
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where m and a s are bosonic mass and s-scattering length, which accounts for the repulsive interactions, and ro is an 
arbitrary spatial scale. Accordingly, the total number of bosons in the gas is given by 

f 2 r N 

N ph = J \ipph(r P h)\ dr ph = ^— , (5) 

where N — J \ip(r)\ 2 dr is the norm of the scaled wave function. 

In the form of Eq. ([3]), the GPE neglects the dipole-dipole interactions between the particles. However, these can 
be readily taken into account, using the same mean-field approximation which is used to derive the GPE. Indeed, 
the local density of the dipole moment in the gas (i.e., the polarization of the medium) is P = d |^(r)| 2 , hence the 
additional electric field generated by the polarization, Erf, is determined by the Poisson equation, V • (Ed + 47rP) = 0, 
which yields E^ = — 47rP = — 47rd | (r) | 2 (this solution for E^ is definitely a unique one for the spherically symmetric 
configurations considered below). Finally, the extra term in the GPE induced by the interaction of the local dipole 
with the collective field, E^, created by all other dipoles is 

-(d-E d )^ = 4^ 2 |^| 2 V. (6) 

Obviously, this term, if added to Eq. ([3|), may be absorbed into a redefinition of the effective scattering length 
accounting for the repulsion between the particles. In the underlying physical units, this amounts to 

a s -> (a s ) off ee a s + md 2 /h 2 , (7) 

where m is the mass of the dipolar molecule. 

Unlike the quantization, the inclusion of the nonlinearity does not break the underlying scaling invariance of the 
Hamiltonian, cf. Eq. (TTJ). Indeed, under the combined transformation, 

r ->• av. ip -> a~ lr tp, fl -> a~ 2 £l. (8) 
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the total energy of the condensate described by GPE (J3]) , 

1 



E D = - 



|V^r-(C/o^ 2 -^ 2 r 2 )|Vr + |V 4 | dr, (9) 



features scaling E — > E/a A ~ D , where D is the dimension (3 or 2). 

It is relevant to mention that a quantum anomaly was very recently predicted in a model described by the GPE 
in the 2D space, for a harmonically trapped gas of bosons interacting through the two-dimensional repulsive delta- 
functional potential The anomaly breaks the specific scaling invariance of this gas, which holds in the mean-field 
approximation [l6j . 

The rest of the paper is organized as follows. In Section II, we briefly recapitulate the description of the 3D and 
2D collapses in the framework of quantum mechanics, extending it through the inclusion of self-similar nonstationary 
solutions, in addition to the known stationary ones. Section III reports the basic results obtained in the 3D case, 
which demonstrate the creation of the previously missing GS by the self-repulsive cubic nonlinearity at Uq > 1/4, as 
well as making the GS normalizable at Uq < 1/4. It is demonstrated too that, in latter case, the inclusion of the 
harmonic trapping potential gives rise to a tristability of bound states. The results obtained in the 2D model with 
the quintic repulsive term are reported in Section IV. In that case, the GS also replaces the quantum-collapse regime 
for Uq > 0. In the case of < — Uq < 1/4 (i.e., in the case of the weakly repulsive potential), the quintic nonlinearity 
also gives rise to the GS with a normalizable wave function. In the latter case, three confined modes are found in the 
presence of the harmonic trap, like in the 3D setting, but the middle mode is unstable, spontaneously transforming 
itself into a breather. The paper is concluded by Section V. 

II. THE COLLAPSE IN QUANTUM MECHANICS 

In the framework of Eq. ([2]), 3D and 2D stationary states, with angular quantum numbers / and m, are looked for 

as 

■03D = exp(-int)Yi m (9,<p)(p(r), (10) 
^2D = exp(—ifj,t + iltp)<f)(r), (11) 

where 9 and ip are the angular coordinates and Y\ m {9, if) is the spherical harmonic. Then, an exact 3D wave function 
in the form of Eq. (JTDJ) can be readily found, as a solution to Eq. @, for Ui = Uq — 1(1 + 1) < 1/4, cf. Ref. 

4>(r) = for- a± cxp (~nr 2 /2) , (12) 



A* = nQ- ff± J, a± = ^±yj\-U h (13) 

the GS corresponding to the smaller value of fi, i.e., <r + . This solution is also relevant for Ui < (the repulsive 
potential), with the respective norm, 

/>oo 

N 3D =4ir cj) 2 (r)r 2 dr, (14) 
Jo 

converging (for <j + ) if Ui > —3/4. 

In two dimensions, the GS solution to Eq. @ exist only for Ui = Uq — I 2 < 0, in the exact form given by Eqs. pT|) 
and |[T2l. but with 



H = n{l-o±),<T±=±y/-U u (15) 
cf. Eq. (fTB")) . Like in the 3D case, the GS corresponds to er + , the 2D norm of this solution, 

N 2U = 2tt / <f> 2 (r)rdr, (16) 



converging if Ui > — 1. 

Past the critical point, i.e., for Ui > 1/4 in 3D, and for any Ui > in 2D, the asymptotic form of solutions for (j>(r) 
at r -> is 0] 

0(r) « 0„ I r " /2 C0S (v^V41° (r/nO) ^ = 3, 
[ cos (y/Ui In (r/r )), D = 2, 
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FIG. 1: Three radial profiles of |x( r >^)l = |V , ( r )l 7 * 1/ ' 2 at t = 0, 0.005 and 0.1 (dotted, dashed, and solid curves), as produced by 
the numerical solution of the radial version of Eq. ((2| in 3D with Uq = 0.27 and fi 2 = 0.1. Note that this value of Uo slightly 
exceeds the critical one, ((To)'^ ' = 1/4- 



with arbitrary constants 4>o and ro (these stationary solutions are exact but unnormalizable ones for ji = CI = 0). 
The infinite number of radial nodes (zeros) in solutions p7|) implies the nonexistence of the GS in these cases, and 
represents the phenomenon of the quantum anomaly. Accordingly, the nonstationary wave function obeying Eq. 
@ is expected to collapse, eventually. The latter conjecture was confirmed by direct simulations. As an example, 
in Fig. 1 we display results of the simulations of the 3D spherically symmetric solution, with I = m = 0. The 
simulations were performed by means of the split-step Fourier method using 2 16 modes. As the initial condition, we 
took tjj(r) — r" 1 / 2 exp(— Sir 2 / 2), which is the exact stationary wave function for Uq = 1/4, i.e., the one at the critical 
point, see Eqs. (|12p and flT3"l) . The snap-shot profiles of r 1 / 2 |^(r)| shown in Fig. [T] testify to the rapid growth of the 
amplitude of the solution at r = 0, which is a signature of the development of the collapse. The growth eventually 
ceases in the simulations, due to the finite mesh size of the numerical scheme. 

A more general version of asymptotic solution (|17|) can be obtained by means of substitution 



which makes Eq. ^ for the 3D time-dependent wave function equivalent to its 2D counterpart: 

i Xt = - (1/2) [d 2 + r- Y d r + (U t - 1/4) r- 2 - fiV] x . 



(18) 



(19) 



Then, omitting the term ~ f2 2 , time-dependent solutions to Eq. (|19p may be looked for in a self-similar form, 
X{r,t) =X(A), where 



A^ln((p 2 /r 2 )(t/r )) 



(20) 



with arbitrary radial and temporal scales, po and To- This substitution transforms Eq. (|19[) with O = into an 
ordinary differential equation, 



dA 



to 
2p§ 



Explicit solutions to Eq. (|21[) can be found in an asymptotic form, corresponding to A — > oo (i.e., r 
t —> oo), which is a spatiotemporal generalization of the stationary solution (JTTJ) : 



(21) 
and/or 



X = 0o cos((l/2)VC/ z -l/4A) 



^Ui - 1/4 r 2 
4Ui + 3 t 



sm 



(l/2)y/Ui - 1/4A 



22) 



+ V/C/Z-1/4COS ({l/2y/Ui - 1/4A) 
the second term being a small correction to the first one. Through substitution (|18l) . solution (l22l pertains to the 3D 



case, and, on the other hand, expression (|2"2l with ^JUi — 1/4 replaced by \fU~i directly applies to the 2D situation. 
With regard to Eq. (|20|) . the meaning of solution (|22p in both the 3D and 2D cases may be understood by noting 

that the number of nodes of the wave function [i.e., the number of the corresponding zeros of cos ^(1/2) y/Ui — 1/4A^] 

which pass a point with fixed r grows with time as n w \fU~i ~ V4/ {2n) [in (r^/r 2 ) + In (t/to)] . This actually 
implies the development of the quantum-anomaly structure in time, which expands outwards. 
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III. THE 3D GROUND STATE CREATED BY THE CUBIC SELF-REPULSIVE NONLINEARITY 
A. The ground state in the absence of the external trap 

Spherically symmetric stationary solutions to the GPE in the form of Eq. ([3]) are looked for as 

^(r,t) = e-^V- 1 X3D(r), (23) 

with real function \ obeying equation 

^X3D =~\ [X3D + {Uor- 2 - ttV) X3D ] + r- 2 X 3 3D . (24) 



The expansion of solutions to Eq. (|24[) at r — > is 

XsdW = + Xir s/2 , 8=1 + Vl + 8C/o, (25) 

where \i is a f ree constant. For any Uq > 0, the cubic nonlinearity supports finite-norm states without the help of 
the external trapping potential (O = 0). Indeed, in this case the asymptotic form of the solution with [i < at r — > 00 
is X3D = Xo ex P ( — V - 2/ir) , with some constant \o- Combining it with the asymptotic form (|25[) valid at r — > 0, one 
may use, as the simplest analytical approximation, the following interpolation for the 3D modes which represent the 
GS with given norm N: 

«-' = # e -"V-^., = -i(^) 2 . (26) 

The small term X1W 2 from Eq. ([25| is not included into this approximation, and the norm was calculated by the 
substitution of analytical expression (|2"o| into integral (fT4")l . Actually, in the limit of /1 — > —0, Eq. (121)1) gives an 
asymptotically exact solution (rather than being simply an interpolation), and its limit form corresponding to /i = 
and N = 00 is an exact (although unnormalizable) solution. 

Numerical solutions of Eq. (|24l) were found by means of the shooting method. A typical example of the GS solution, 
along with the respective approximation (f2l)|). is displayed in Fig. [U[a) for Uq = 0.8, which is essentially larger than 
the critical attraction strength, (Z7o)^ D ' ) = 1/4, beyond which the GS does not exist in the framework of the linear 
Schrodinger equation @. Further, Figs. HJb) and [He) represent the entire family of the solutions for two values, 
Uo = 0.8 and 0.1, which are, respectively, larger and smaller than 1/4. The conclusion is that, in the nonlinear model, 
the GS exists for all values of Uq and N. Thus, the self-repulsive cubic term completely suppresses the quantum 
collapse in the 3D space and offers an alternative solution to the quantum-anomaly problem by creating the GS 
where it does not exist in the linear Schrodinger equation. 

The analytical approximation d26l) suggests an estimate for the radial size of the GS created by the repulsive 
nonlinearity: 



R (3D) _ 4tt 
GS "iV3D 



'0 



It is relevant to rewrite this estimate in terms of the physical units, as per Eqs. (j4j, <j5j> , and ([7]): 



(«<??>), sn n a3 = 2 J^±pfl^, ( 28) 

V / P h m{U ) ph 

which gives the radius of the GS in terms of the physical parameters of the model. Note that arbitrary spatial scale 
ro, which was used in rescalings (j4]) and ([5]), does not appear in Eq. ([5]) (actually, it cancels out). It is natural that 
Rqs shrinks to zero at iV — > 0, which implies the transition to the collapse in the framework of the linear Schrodinger 
equation. 

The existence of the GS with the finite norm at Uq < 1/4 in the nonlinear model with £1 = is worthy noting 
too, as the corresponding stationary solutions to the linear equation, see Eqs. (|10[) and (|12j) . are unnormalizable 
for O = 0. Lastly, simulations of the radial version of Eq. with arbitrary random perturbations added to the 
stationary solutions (not shown here) demonstrate that the entire GS family is stable. Those simulations did not test 
the stability of the GS against tesseral (spherical-symmetry-breaking) perturbations, but the repulsive character of 
both the contact and dipole-dipole interactions (in the mean-field approximation) makes the presence of a symmetry- 
breaking instability implausible. 
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(b) (c) 




FIG. 2: (a) A typical example of the ground state, shown in terms of x( r ) = r IV'MIj hi the 3D nonlinear model without the 
external trap (Q — 0), for Uo = 0.8 and /i = —0.225. Panels (b) and (c) display curves fi(N) for the ground-state families 
with Uo = 0.8 and 0.1, respectively. In all the panels, the solid and dashed curves depict, severally, the numerical results and 
analytical approximation (|26|l . In particular, the latter predicts N(fi = —0.225) = 5.30 for Uo = 0.8 [the case shown in (a)], 
while the numerically found counterpart of this value is iV num (/i = —0.225) = 6.26. The convergence of the numerical and 
analytical curves for iV(/z) at fj, — oo corresponds to the fact that Eq. (|26p gives an asymptotically exact solution in this 
limit. 

It is relevant to mention that, as follows from Eq. (|25|) . the (scaled) energy of the GS, calculated, as per definition 
(IH]), in a regularized form, i.e., for r > g — > 0, contains a diverging term, 

E 3B = irU (2 - Uo) (2Q)- 1 , (29) 

which may be removed by means of the renormalization, cf. Refs. [H0|- The vanishing of E^d at Uo = 2 seems to be 
a formal peculiarity, rather than a physical feature of the model. 

B. Effects of the harmonic trap: the tristability 

The addition of the harmonic trap deforms the GS family reported above for fi = 0. As shown in Fig. [31 at 
Uo < 1/4 the nonlinear model with ft > supports two additional families of 3D confined modes, thus featuring a 
tristability. The lowest curve in Fig. GH^a) represents the deformed GS branch produced by the numerical solution 
of Eq. (|24[) . while two upper branches represent a nonlinear deformation of the exact solutions generated by linear 
equation © in the form of Eqs. (fTU)) . (IT21 . and (fT3")l (note that in the linear limit, i.e., at N — > 0, the GS created 
by the nonlinearity disappears by falling to \i —> — 00). The upper branches merge and disappear at Uq —> 1/4. The 
stability of all the three families was verified by direct simulations of the radial version of Eq. §3§ (not shown here). 
In the limit of large N, all the three branches in Fig. [3Ja) asymptotically approach an expression predicted by the 
Thomas-Fermi approximation, 

(N TF ) 3D » (l6\/27r/15fi 3 ) a* 5/2 - (30) 



C. 3D vortical modes 

Nonlinear solutions for 3D states carrying the angular momentum may be considered by means of the averaging in 
the angular variables, \Yi m (9, tp)\ 2 Y lm (0, ip) -> c im Y" (m (9, tp) with q m = sin 9d9 J^ n \Y lm (9, <^)| 4 dip (in particular, 
c io = 9/5, C\i — 6/5). In this approximation, ansatz 

^ = ^ /2 Yim (9, tp) e-**r- x x(r), (31) 

cf. Eq. (|2U)l , leads to Eq. ([M]) with x replaced by x, and Uo substituted by Uo — I {I + 1)- In Fig. 0] we compare 
numerically found curves n(N) for the modes with (/, m) — (0, 0) and (1, 0), fixing Uo — 2.2 and Q = 0. In this figure, 
the mode corresponding to (l,m) = (1,0) was generated using the spherically symmetric equation, with the norm 
modified as per the averaging approximation, see Eq. (I31[) . 
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FIG. 3: (a) Curves fi(N) for three families of the 3D confined modes at fi 2 =0.1 and Uo = 0.2. (b) Profiles of typical modes 
belonging to the upper, middle, and lower branches in (a) are shown by the solid, dashed, and dotted curves, respectively. In 
panel (b), x = r \ % l } ( r )\ f° r the solution belonging to the lower branch in (a), while \ = r a± \ip(r)\ for the two upper ones, with 
the same a± as in Eq. (I13|l . The profiles are shown in (b) for N = 10.2, the respective values of \x being 0.517, 0.503, and 
0.460. 




FIG. 4: The /J.(N) curves at Q 2 — and Uo = 2.2 for 3D modes with angular quantum numbers (l,m) = (0,0) and (1,0), 
obtained by means of the approximate ansatz (|31[l , 

IV. 2D GROUND STATES CREATED BY THE QUINTIC SELF-REPULSIVE NONLINEARITY 
A. The ground state in the absence of the external trap 



As said above, the GPE in the form of Eq. (J3|) is relevant in the 2D case too. However, the 2D norm of the solution 
with asymptotic form ~ at r — > 0, which follows from this equation [see Eq. ([26]) ]. diverges. In other words, the 
cubic self-repulsion is not strong enough to prevent the collapse in the 2D setting. On the other hand, the GPE may 
also include the quintic repulsive term accounting for three-body collisions, provided that the collisions do not give 
rise to conspicuous losses [17] . The 2D axisymmetric GPE with the dominating quintic term is 



(1/2) (Vw + r~ Vr + Uir~ 2 - Q 2 r 2 ) i> + |^| 4 if>, 



(32) 



[recall V~i = U — I 2 , if vorticity I is present, see Eq. (fTTj) ; unlike the 3D case, the vorticity may be kept in the 2D 
analysis based on radial equations]. In Eq. (j32j), we neglect inessential cubic terms, including the one which may be 
generated by the dipolc-dipole interactions, as in Eq. (J7|). If the cubic terms are kept, they do not significantly affect 
the results presented below [in particular, they do not alter the first term in expansion (f36f at r — > 0]. 
The total energy of the BEC described by 2D equation (|32|) is 



E2D — 77 



|VV>| 2 - (U Q r- 



dr, 



(33) 



cf. Eq. ([5]). This energy features the invariance with respect to the scaling transformation, E-2D — > -E^d/c*, with the 
difference from the 3D case in that the wave function is transformed as per ip — > a -1 / 2 ?/;, cf. Eq. (|SJ). 
Stationary solutions to Eq. (13^1) are looked for as 



V-M) 



X2i)(r), 



(34) 
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cf. Eq. ([23]) . which yields an equation for x(r): 

1 

MX2D = -7; 



X2D +[Ui + l) r-\ 



+ r~ 2 X 5 , 



where we set f2 = 0, cf. Eq. (|24p in 3D. The expansion of the solution to Eq. ([33)) at r — > is 



1/4 



(35) 



(36) 



where s = (1/2) (l + y/5 + 16{7/), and Xi is an arbitrary constant, cf. Eq. (|23|) in the 3D case. The solution with a 
finite norm exists at Ui > —1/4, representing, at Ui > 0, the suppression of the collapse and creation of the GS by 
the quintic nonlinearity. 

As concerns the vorticity, in the usual situation the amplitude of the corresponding mode must vanish at point r = 0, 
where the phase cannot be defined. However, in the present case the solution features a different solution to the phase- 
singularity problem: instead of vanishing, the amplitude diverges at r — > - as \ip(r)\ « [(1/2) (Ui + 1/4)] 1 / 4 r^ 1 / 2 , 
according to Eqs. (|34[) and (|36p. This possibility is akin to the fact that the Bessel equation for amplitudes of vortex 
modes with topological charge I gives rise to two solutions, J±\i\ (const • r), the one with index —\l\ being usually 
omitted as a singular one. However, in the present case the singular solution is acceptable, as it provides for the 
convergence of 2D norm (|16p . 

The energy of the GS, if calculated at r > g —> 0, contains a divergent term, E2D = T^yUi + 1/4 (1 — 2U{) (3-\/2f?) , 
cf. the similar term (1291) in the 3D setting. This term, which may be removed by the renormalization procedure, and 
its vanishing at Ui = 1/2, do not seem to be physically significant features of the model. 

Combining the 2D asymptotic form (|36j) . valid at r — > 0, and its counterpart at r — > 00, \2D ~ Xo ex P ( — V — 2/ir) , 
and making use of definition (|16[) for the 2D norm, we obtain an analytical interpolation formula for the GS family, 
in the absence of the external trap (O = 0): 



^2D 



Ui 



-,1/4 



e-*#r-V 2 e- 



(37) 



cf. approximation ([26]) in the 3D case. Similar to the situation in the 3D case, Eq. ([37]) gives an asymptotically 
exact solution (rather than a mere interpolation) for /1 — > —0, and an exact solution with the infinite norm at /i = 0. 
The approximation (|3T[) makes it possible to define the radial size of the two-dimensional GS created by the quintic 
nonlinearity, cf. Eq. ([27| in the 3D case: 



R, 



(2D) 
GS 



2tt 



T2D 



0) 



(r) 



r dr 



2D 



^2(^0 + 1/4) 



(38) 



Note that the quintic term supports the GS in 2D even at < — Ui < 1/4, when the central potential is repulsive. 
The correctness of this counter-intuitive conclusion is corroborated by the above-mentioned fact that the analytical 
approximation (|37p gives the asymptotically exact solution for fi — ¥ 0, including the case of < — Ui < 1/4. 

An example of the stable GS, and curves n(N) for the GS families in 2D are displayed, along with the analytical 
approximation (13"T1) . in Fig. [S] The fi(N) curves are shown for both signs of the central potential (and 1 = 0), 
Uq = —0.18 and Uo = 0.05. Simulations of perturbed solutions within the framework of the radial version of Eq. 
([32]) confirm the stability of the GS families (not shown here) . Although these simulations do not include azimuthal 
perturbations, the repulsive sign of the nonlinearity makes it evident that these perturbations will not give rise to an 
instability 



B. Effects of the harmonic trap 



As shown in Fig. EJa), in the presence of the external trap (Q > 0), three different stationary solutions can be 
found at —1/4 < Ui < 0, when the corresponding 2D linear Schrodinger equation gives rise to two exact solutions, in 
the form given by Eqs. ([TTj) . (TT2"j) . and flTSl) . As well as in the similar situation for the 3D case, which is displayed in 
Fig. Eta), two upper branches in Fig. [S] represent nonlinear deformations of the exact linear solutions, with values of 
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(a) (b) 




FIG. 5: (a) The radial profile of the ground state in the 2D model with the quintic nonlinearity and £1 = 0, for Uo ~ 0.05 and 
H = —0.1867. (b) Curves n(N) for the GS families with Uo = —0.18 and Uo = 0.05. In both panels, the numerical results 
and the respective analytical approximation (I37[) are depicted by the continuous and dashed curves. The convergence of the 
numerical and analytical curves for N(fi) at /i — s> —0 corresponds to the fact that Eq. (|37[) gives an asymptotically exact 
solution in this limit. 



(a) (b) 




FIG. 6: (a) Curves fl(N) for three stationary modes found in the 2D model with the quintic nonlinearity for Uo = —0.18 and 
Q, 2 = 0.10, cf. the similar picture for the 3D cubic model displayed in Fig. (b) An example of the breather generated by 
the solution with N = 0.3453 belonging to the unstable middle branch from panel (b). The plot shows A = \\(r = 0,i)| 2 as a 
function of time, with \ defined as X2D in Eq. ()34f) . 

fi approaching those given by Eq. (1151) in the limit of ./V — > 0. The lowest branch represents the GS corresponding 
to approximation (f3"T|). which is additionally deformed by the trapping potential. In the limit of \l — > oo, both the 
bottom and top branches in Fig. [S{a) asymptotically approach the Thomas-Fermi limit, (-/Vtf) 2 d ~ (47r/3f2 2 ) /i 3 / 2 , 
cf. Eq. P0]h 

In spite of the overall similarity to the 3D case, it is worthy to note that the shape of the middle branch in Fig. 
[f^a) is drastically different from its counterpart in the 3D case, cf. [3Ja). Moreover, direct simulations of Eq. (j3"2")l 
demonstrate that, while the solutions corresponding to the top and bottom branches in Fig. IHIa) are stable (not 
shown here in detail), the middle-branch solutions are not. In the simulations, they are spontaneously transformed 
into robust breathers featuring long-period oscillations, see an example in Fig. [6f b) . 

V. CONCLUSION 

We have demonstrated that the mean-field repulsive nonlinearity suppresses the 3D quantum collapse induced by 
the central attractive potential, — (Uq/2) r~ 2 , which can be realized in the ultracold gas of dipolar molecules attracted 
by the central charge. The dipole-dipole interactions were also taken into account, resulting (in the framework of the 
mean-field approximation) in a redefinition of the scattering length which accounts for the contact repulsion. The 
nonlinearity creates the GS (ground state) in place of the collapse regime. For Uq < 1/4, when the respective 3D 
Schrodinger equation does not lead to the collapse, the inclusion of the harmonic trap gives rise to the tristability. 
The cubic repulsion is not strong enough to prevent the 2D collapse, but the quintic term is sufficient for this purpose. 
It also gives rise to the GS which replaces the quantum collapse in the 2D space, and to similar modes carrying the 
angular momentum (which feature the amplitude diverging, rather than vanishing, at r — > 0, while the total norm of 
the vortical mode converges). A difference from the 3D case is that, in the presence of the harmonic trap, one of the 
three confined modes supported by the weakly repulsive central potential with < — Uo < 1/4, in the combination 
withe quintic nonlinearity, is unstable, transforming itself into a breather. A counter-intuitive finding is that the 2D 
self-trapped mode exists even in the case of the weakly repulsive potential, while the harmonic trap is absent. 



10 



This work suggests continuations in several directions. If the orientation of the dipoles in the 3D space is fixed 
by an external uniform field, the central charge induces an axisymmetric potential, U = — (Uq/2) r~ 2 cos 8, and it 
may be interesting to study the possibility of the replacement of the corresponding anisotropic collapse (cf. Ref. [7|) 
by a GS. It is relevant to mention that the dipole-dipole interactions between the bosons may give rise to a specific 
mode of the nonlinear collapse in BEC 19], in the absence of the contact repulsion, which suggests to consider an 
interplay of this nonlinear mode with the linear collapse (although the dynamics of the collapse is not described by 
the mean- field approximation [l2j). Another interesting extension may be the study of higher-order nonlinear states, 
such as 3D vortical modes. 

A challenging problem is to extend the analysis to fermion gases. In that connection, it is relevant to mention 
that the fall-onto-the-center effect was studied for interacting fermion pairs, in the framework of the Bethe-Salpeter 
equation [20l |. 
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